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Introduction 

Most geometric modeling systems use either polynomial or rational 
functions to represent geometry. In such systems most comput- 
ational problems can be formulated as systems of polynomials in 
one or more variables. Classical elimination theory can be used 
to solve such systems ( Refs. 1,2, 3, 4, 5 ). In this paper we 
summarize Cayley's method of elimination and show how it can best 
be used to solve the curve/curve intersection problem. 
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Res(P,Q) is a polynomial in x i> 
three conditions are equivalent 

1. Res (P,Q) = 0 

2. P and Q have a common root 

3. The mxm system of linear homogeneous equations 
given by r 1 r x 0-J = o has a solution. 
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Furthermore , there ( usually ) exist polynomials 
F (xj , — ’X n _j ) , G(x x ,. . . ,x ) such that if Res(P,Q) 
= 0 for some fixed x x> . .. » x n i then the common root 


of P and Q is given by 


= -F/G. The functions F 


and G are derived by Cramer's rule from the system 

[c. .] rxV = o • 
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Example. 
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Clearly detfc^j] = 0 , and from 
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it follows that -4+4x = 0 and thus x = 1. 
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Curve/Curve Intersection Using Elimination 

Let C t (s ) = (Xj (s)fYi ( s ) ) and C 2 (t) = (x 2 (t) ,y 2 ( t) ) be two 
rational parametric curves , both defined on [0,1]. There are two 
ways to use elimination to find the intersection points of C x and 
C 2 • 

Method 1. Implicitize one of the curves. 

Steps : 

a. Use elimination to convert C x (s) to its implicit represent- 
ation P(x,y) = 0. In the process we obtain s = -F(x,y)/G(x,y) 
This is an interesting problem in its own right. 

b. Substitute x 2 (t) , y 2 (t) into P(x,y) =0 to get Q(t) = 0 . 

c. Find the roots t-j of Q(t) within the range [0,1]. 

d. For each t-j , use x 2 (t) , y 2 (t) to obtain the corresponding 
point (V y i) on C 2 . 

e. Use (x i ,y 1 - ) and s = -F/G to get s i . If s i is in [0,1] , then 
s i , t. , (x.. ,y n - ) is a solution ( intersection point ). 

Example. 

Let Cits) = ( (1-s 2 )/(l+s 2 ) ,2s/(l+s 2 ) ) and C 2 (t) = (t,t 2 ). C x is a 
circular , C 2 a parabolic arc , both in the first quadrant. 

a. From x= ( 1-s 2 ) /(1+s 2 ) and y=2s/(l+s 2 ) , it follows that 
f(x,y,s) = (x-1) +(x+l) s 2 = 0 

g(x,y,s) = y-2s+ys 2 = 0 . Eliminating s , we obtain : 

-2x+2 -2y 

[ Cij ] = 

|_ -2y 2x+2_ 

Expanding det^ c ij^ and setting it equal to zero yields P(x,y) 

= x 2 +y 2 -l = 0. And from (-2x+2)-2ys = 0 , we get s = (l-x)/y. 

b. Substituting (t,t 2 ) into P(x,y) yields Q(t) = t 4 +t 2 -l = 0. 

c. The only root in [0,1] is (approx.) t = 0.786 . 

d. From C 2 (t) , we have x = 0.786 , y = 0.618 . 

e. From s = (l-x)/y we obtain s = 0.346 . 

Method 2. Subtract the coordinate functions. 

Steps : 

a. Form P(s,t) = Xj(s)-x 2 (t) = 0 and Q(s,t) = y } (s)-y 2 (t) = 0. 
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b. Eliminate one of the variables , say s , to get Res(P,Q) = 
R (t ) = 0 and s = -F(t)/G(t). 

c. Find the roots t.. of R(t) in the range [0,1] . 

d. For each , use s = -F/G to get . if is also in 
[0,1] , use (x 2 (t) ,y 2 (t) ) to get ( x y ) . 

Example. 

Same as above. 

a. 0 = P(s,t) = ( 1— s 2 ) / ( 1+s 2 ) - t = (l-t)-(l+t)s 2 
0 = Q(s,t) = 2s/ (1+s 2 ) - t 2 = -t 2 +2s-t 2 s 2 


b. 

R(t) 


2 ( 1-t) -2t 2 

-2t 2 2 ( 1+t) 


= t 4 +t 2 -1 = 0. 


And 2(l-t)-2t 2 s = 0 implies s = (l-t)/t 2 . 

c. Solving for t yields (approx.) t = 0.786 

d. Substituting into s = (l-t)/t 2 yields s = 0.346 


Comparison of the Two Methods 

To our knowledge only method 1 has been mentioned in the CAD/CAM 
literature ( Refs. 2,3 ). But method 2 is a more straight- 
forward approach. Furthermore , it is computationally simpler , 
since the elements of the Cayley matrix are one variable instead 
of two variable polynomials. We implemented and tested both 
methods and found method 2 to be more efficient. We used six 
pairs of curves , representing mixtures of lines , circles , and 
cubic arcs. Several examples had multiple intersection points. 

For all six cases method 2 required less CPU time than method 1. 

The average time ratio of method 1 to method 2 was 3.13:1 , the 
least difference was 2.33:1 , and the most dramatic was 6.25:1 . 

Conclusion 

Both of the above methods can be extended to solve the 
surface/surface intersection problem. That is the direction of 
our current research. 
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